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ABSTRACT 

High resolution observations of Sgr A* have revealed a wide variety of phenomena, ranging from intense 
rapid flares to quasi-periodic oscillations, making this object an ideal system to study the properties of low 
luminosity accreting black holes. In this paper, we use a pseudo-spectral algorithm to construct and evolve a 
three-dimensional magnetohydrodynamic model of the accretion disk in Sgr A*. Assuming a hybrid thermal- 
nonthermal emission scheme, we show that the MHD turbulence can by itself only produce factor of two 
fluctuations in luminosity. These amplitudes in variation cannot explain the magnitude of flares observed in this 
system. However, we also demonstrate that density perturbations in the disk do produce outbursts qualitatively 
similar to those observed by XMM-Newton in X-rays and ground-based facilities in the near infrared. Quasi- 
periodic oscillations emerge naturally in the simulated lightcurves. We attribute these to non-axisymmetric 
density perturbations that emerge as the disk evolves back toward its quiescent state. 
Subject headings: accretion — black hole physics — Galaxy: center — instabilities — MHD — relativity 
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1. INTRODUCTION 

Accreting black holes of all masses, from stellar-mass sys- 
tems in X-ray binaries to supermassive objects in active galac- 
tic nuclei (AGN), are highly variable, exhibiting a wide vari- 
ety of outbursts from simple flares to relativisti c jets. Peri- 
odic o r quasi-periodic signals in these outbursts (Ivan der KlisI 
l2006l) . if connected to the orbital period at the black hole's 
innermost stable circular orbit (ISCO, also known as the 
marginally stable orbit), provide ideal probes of the spacetime 
in that region (see Melia et al. 2001a; Psaltis 2003). How- 
ever, using observations of the emission from the ISCO as 
probes requires understanding the dynamics of accretion near 
the black hole and the origin of these oscillations. 

Although these quasi-periodic oscillations (QPOs) are well- 
established observationally, particularly in stellar-mass sys- 
tems, their origin is still unclear. Quasi-periodic variabil- 
ity has emerged from theoretical work only under restric- 
tive a ssumptions, either analytically in idealized disks (e.g., 
lKatoll200U and references cited therein), or numerically in 
hydrodynamic (e.g., Milsom & Taam 1997) and magnetohy- 
drodynamic (MHD; e.g.. iTagger & Mehall2006h simulations 
with large-scale magnetic fields. However, global calcula- 
tions including the effects of the magnetorotational insta- 
bility (MRJ) and MRJ-driven turbulence have not yet pro- 
duced quasi-periodic signals^ (e.g., Hawley & Krolik 200T|; 
lArmitage & Revnoldsl2003l: lArras et alj|2006, and references 
cited therein). In addition, it is possible that the origin of these 
oscillations is different for different accreting black hole sys- 
tems. 

At the Galactic center, the compact radio source Sagitar- 
ius A* (Sgr A*) derives its power from accretion onto a 
~ 3.4 X 10^ Mq black hole (ISchodel et allllool iGhez et alJ 
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l2004l) . Occasionally, rapid flares are observed from this object 
in the X- rays (Baganoff et a l. 2001) and in the near infrared 
(NIR; see lGenzei et alj2003l) . During these flares, there is ev- 
idence for a quasi-periodic modulation of the X-ray emission, 
with a period of ^ 17 mins that is comparable to the orbital 
period at the ISCO of a Schwarzschild black hole with the 
mass of Sgr A*. 

The spectra and polarization of the source during the qui- 
escent state of Sgr A* suggest that the long-wavelength emis- 
sion is due to synchrotron radiation, whereas the more ener- 
getic X-ray emission is due to inverse Compton scattering and 
thermal bremsstrahlung (see, e.g., Melia 1992; Naravan et aQ 
Il995h . The rapid increase in the NIR and X-ray fluxes during 
the flares can then be accounted for by a transient acceler- 
ation of the energetic electrons resp onsible for the emission 
dMarkoffet al.ll2001l; lYuan et alj |2004) or by a rapid mcrease 
of the accretion rate. Such an increase can be caused either 
by the highly variable nature of the accretion flow or by an 
external surge of matter In this paper, we will address the 
possibility that the X-ray flares are due to a rapid increase of 
the accretion rate, including situations in which this is induced 
by infalling clumps of plasma. 

We will us e the pseudo-spectral algorithm of Chan et alj 
(I2005L l2006allbl) to simulate the effects of the magnetohydro- 
dynamic turbulence on the accretion disk that surrounds the 
black hole. We will focus on the long-wavelength emission 
of Sgr A* and, therefo re, consider only synchrotron emission 
from a hybrid plasma (lOzel et al.ll20 00). For a study of spec- 
tral properties including synchrotron emission/absorption, 
free-free emission/abs orption, and Compt on/inverse Compton 
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scattering, we refer to'Ohsuga et aTl ( 1200 5). We calibrate our 
simulations using the quiescent spectrum of Sgr A* and aim 
to simulate the dynamical evolution of the system caused ei- 
ther by the MHD turb ulence or by an exter nal surge of mat- 
ter. In agreement with lGoldston et aP (l2005h . we find that the 
MRI-driven turbulence alone cannot produce variations in the 
radiation flux that are as large as the observed flares. We also 
do not find any significant quasi-periodic oscillations during 
the quiescent-state flux. 
To produce the observed flares, we locally perturb the 
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density of the disk to simulate the effects of "clumpy 
materi al" raining do wn o nto it from the larg e-scale 
flow (Falcke & Melid [19971 ICoker. Meha. & Falckd [T999t 
llagger & Melia 2006ft. Here, again, observations help to con- 
strain the perturbation, particularly its location in the disk. 
Most of Sgr A*'s luminosity is emitted at the mm/sub-mm 
spectral excess, suggesting for our calculation th at the accre- 
tion disk extend s out to 5-25 Schwarzschild radii (Melia et"an 
llOOO; Melia & Falckdl200lh . Our density perturbations are 
introduced locally within the disk, rather than from a global 
simulation of the infall. However, even with this limitation, 
we find that we have enough freedom to construct a den- 
sity perturbation that produces flares qualitatively similar to 
those observed in Sgr A*. The simulated flares exhibit quasi- 
periodic oscillations. If such density perturbations are indeed 
the cause of black-hole outbursts, this model may be used 
to not only study the properties of spacetime near Sgr A*, 
but also the near-horizon environment in other low-luminosity 
AGNs. 

In ^of this paper we present our modifications to the hy- 
brid thermal-nonthermal emission model of Ozel et al. (200Q) 
by introducing a cooling break to the nonthermal component 
and a more reasonable treatment of the low energy cut-off in 
the electron distribution. In ^ we focus on the observational 
constraints on both our emission model and the conditions at 
the outer boundary of our simulated disk. In ^ we sum- 
marize the physical setup of our simulations based on obser- 
vations. In ^ we present our results on the quiescent state 
and compare the properties of the magnetohydrodynamic tur- 
bulent plasma with those of previous numerical studies. We 
present the results of simulations involving a density perturba- 
tion in the disk in ^ We also discuss accretion rates, the disk 
morphology, and lightcurves from the quiescent and perturbed 
simulations in this section. We summarize the limitations of 
our simulations in ^ We conclude with a discussion of the 
broader impact of these results in ^ 

2. EMISSION MODEL 

In this section, we first describe how to incorporate special- 
and general-relativistic effects for the radiative transfer in 
a pseudo-Newtonian gravity. We then propose a hybrid 
thermal-nonthermal model for the electron distribution and 
consider its synchrotron emissivity from a region near the 
ISCO. 

2.1. Radiative Transfer in Pseudo-Newtonian Gravity 

We will use cylindrical coordinates (r, (f>,z) throughout this 
paper. Because the vertical size of our computational domain 
is small compared to the disk radius, we will ignore gravity 
in the vertical direction. We denote the observed frequency 
at infinity by i'q. The corresponding frequency measured 
by a stationary observer in the local free-falling frame in a 
Schwarzschild spacetime, i.e., without gravitational redshift, 
is given by the transformation 
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where rs = IGM/c^ is the Schwarzschild radius, and G, M, 
and c are the gravitational constant, mass of the central black 
hole, and speed of light, respectively. The flux density ob- 
served at Earth is 



(2) 



where D w 8.5kpc is the distance to the Galactic center and 
is the specific intensity measured in the local free-falling 
frame. The area element, dA, with general-relativistic correc- 
tion, is 

cos/ 

dA = — rdrdcj) , (3) 

where / is the inclination angle between the disk axis (i.e., the 
z-axis in our simulations) and the line of sight. 

The specific intensity is computed in a frame comoving 
with the plasma. The transformation between the comoving 
frame frequency and the local free-falling frame frequency is 



1 -/?cosa(, 



(4) 



and the corresponding transformation for the specific intensity 
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where cos a[, is the cosine of the angle between the velocity 
V and the line of sight. If v = v/|v| is the directional vector 
of velocity and i = -f sin /sin 0-0 sin /cos cos / points to 
the observer from the disk, then 



cos a„ 
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Note that for azimuthal velocity we use V0 — jfSc to take into 
acc ount special-relativistic e ffects in pseudo-Newtonian grav- 
ity dAbramowicz et aLlll996h . so that 
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is always less than unity. 

For simplicity, we assume a time-independent transfer 
equation 

dli, 

—TT^jiy-aJu (8) 

dl 

for each snapshot of our simulations, where / is the line ele- 
ment along the ray (parallel to the directional vector i), and ji, 
and ai, are the emission and absorption coefficients, respec- 
tively. The source function of hybrid synchrotron emission, 
i.e., the ratio of the emission to the absorption coefficient, 
Su = jv/dv, is not equal to a blackbody at the local tem- 
perature. This is especially true at wavelengths l onger than 
the pe ak of the radio/NIR spectrum of the source dOzel et alj 
I2OOOI) . On the other hand, near the peak of the radio/NIR 
spectrum, the difference between the correct source function 
and the blackbody function is negligible as long as the non- 
ther mal electrons are a small fraction of the thermal popula- 
tion (l Ozel et ani2000l) . Furthermore, at even shorter wave- 
lengths, where the emission is optically thin, the source func- 
tion does not enter the calculation. Because evaluating the 
hybrid source function at every wavelength for our particular 
electron distribution is a time-consuming numerical step, we 
will focus our attention to wavelengths comparable or shorter 
than the peak of the radio/NIR spectrum and approximate the 
source function with the blackbody function. The transfer 
equation, written in terms of physical depth, takes the form 



dlv 
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Fig. 1 . — The electron distribution function ^(7) for different choices of 
The other model parameters are fixed as shown in the legend. At the cooling 
break 71,, the power-law index increases by 1. 

Note that this approximation has only minor effects on the 
optically-thin portion of the radiation spectrum. This equation 
can be integrated numerically along rays parallel to the line of 
sight by the first order forward difference equation 



r«+l 



AI 



(10) 



through the computational domain. The superscript n denotes 
the steps and Al = (Az/cos/) is the finite difference "line 
element". For each ray, if /" exceeds Bi, before integrating 
through the whole domain, we simply stop and set Ii, equal to 
the local blackbody value. 

2.2. Electron Distribution 



lOzel et al] (l2000h proposed a hybrid thermal-nonthermal 
model for synchrotron radiation in low-luminosity AGN. 
When applied to Sgr A*, the hybrid model shows better agree- 
ment with observations compared to a pure thermal model. 
However, recent polarization measurements of Sgr A* show 
that the rotation measure in the emission region is small, 
which suggests that the low-energ y component of ele ctrons 
is relativistic and likely thermal ( Marrone et al.ll2006l) . We 
modify the hybrid model so that the nonthermal distribution 
contributes only at high energies. That is, for some critical 
value 7c, we assume 



A^r-Y^ - / ^*(7) , 7<7c, 
^(^^ - \ A^p,(7) , 7>7c, 



(11) 



where 7 is proportional to the (ultra-relativistic) electron en- 
ergy. When most of the electrons are thermalized, our distri- 
bu tion function is a close approximation to the hybrid model 
of lOzel et al.l (l2000h . Otherwise, our modified hybrid model 
keeps the low-energy electrons in a thermal distribution. 

In order to obtain analytical expressions for the normaliza- 
tion, we take the domain of 7 to be all positive real numbers, 
and use the ultra-relativistic Maxwell-Boltzmann distribution 



Mh(7) =A^c^exp 
7c 



7c-7 



(12) 



where A^c is the normalization and 9^ = ksTf/m^c^ is the di- 
mensionless electron temperature. We use to denote the 
Boltzmann constant, to denote the electron temperature, 
which is assumed to be equal to the plasma temperature T, 



and nie to denote the electron mass. Provided 9^ ^ I, the 
ultra-relativistic approximation will not introduce an appre- 
ciable error. Note that the maximum of the thermal distri- 
bution is located at 7 = 29^. It is convenient to introduce a 
parameter 

C = 7c/2^.>1, (13) 

and specify it in our model. 

We also introduce an additional parameter for the cooling 
break. The broken power law is given by 



A^p,(7)=A^c^ 7b7cV7-'^' , 7>7b 



(14) 



The power-law index s describes the spectrum of injected 
electrons and is believed to be greater than unity. Because the 
synchrotron cooling timescale is proportional to 7^^, the elec- 
trons in the high-energy tail cool more rapidly. The parameter 
7b therefore controls the location of this cooling break. 

The symbol A^c in the above equations denotes the normal- 
ization of the distribution at 7c. Integrating the distribution, it 
is easy to show that 
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+ -7b 
s 



^ ,(15) 



where p is the mass density and nin is the mass of hydrogen. 

Figure [T] shows the electron distribution function A^(7) for 
different choices of ^ with fixed density = p/niH and tem- 
perature. The function approaches a pure thermal distribu- 
tion as — > 00. As we decrease (, more and more electrons 
are placed in the nonthermal tail. Note that the spectral in- 
dex of the power-law component increases by one at the lo- 
cation of the cooling break 7b. This modified hybrid thermal- 
nonthermal model will be used to fit the spectrum of Sgr A* 
up to the X-ray band in ^ 

2.3. Synchrotron Radiation 
The emission coefficient in equation (|9]) is given by 

V3e^ 



-lyBsmaB I N{j)F{x)d^ . 
2c Jo 



(16) 



where i/b = eB/2TimeC is the non-relativistic cyclotron fre- 
quency, and aB is the angle between the magnetic field and 
the line of sight in the comoving frame. Here we use e to de- 
note the electron charge and B to denote the magnitude of the 
magnetic field. The function F{x) is defined by 



(17) 



F{x) = x / K5/j,{z)dz 



(iPacholczykll 1 970h . where x is related to 7 according to x = 
xo/^'^, with 

2iy 

xo = ^ : . (18) 

jt'Bsm aB 

We assume that the sub-grid fluctuation of the magnetic field 
is small. Writing B = B/|B|, we simply use 

cosa^ = B-i (19) 

to obtain the angle between the magnetic field and the line of 
sight. The comoving angle cosine is then given by the trans- 
formation 

cos Qfg - /? 
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Fig. 2. — Dependence on the model parameters of the flux density F,^ observed at Earth. The flux densities are computed from our emission model assuming 
a simplified disk. The top panels [from (a) to (c)] present the dependence on different parameters in the emission model, while the bottom panels present the 
dependence on the parameters of the disk itself. The units for T are K; the units for n,, are cm"'' . The non-monotonic dependence of the spectrum on cos i is 
caused by the combined effects of Doppler boosting, projection, and the pitch-angle dependence of the synchrotron radiation. 



Using xm = the emission coefficient can be written 



V3c 
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(21) 



The partial thermal synchrotron function, (jcm), is given by 
the integral 

r2C 



,2 



dz , 



(22) 



1 

Xm Jq 

where z — j/Oe- It does not have any known analytical form. 
We approximate it by a piecewise power law, i.e., for each 
fixed (, we pre-compute the function numerically as a lookup 
table and carry out linear interpolations on a log-log scale. 
Similarly, the function Js.-yf,{xo,jc) is given by 



-/.v.7b(jC0,7c)^ 



1 



7b 



'''-^'^^^F{x)dx , 



(23) 



where = x^j^^ and Xb — xo/7b^- The integrals are related 
to a class of hypergeometric functions. However, because 
there is no convergent algorithm to compute this specific class 
of hypergeometric functions, we will simply perform the inte- 
gral numerically, for each s and 7b. We use a piecewise power 
law to approximate the first integral as a function of Xg. 

Figure |2] demonstrates how the spectrum depends on the 
different parameters. To highlight the characteristics of the 



emission model, we assume for illustrative purpose that the 
disk is uniform and rotates with a constant (dimensionless) 
azimuthal velocity = 2/3. We also assume an azimuthal 
magnetic field, whose energy density is 10% of the internal 
energy of the gas, and use a disk volume of \tx x lO^^cm^, 
which is of the same order as our computational domain (see 
® . In the first row of panels, we vary different parameters of 
the electron distribution. It is interesting to note that in panel 
(a), the nonthermal tail depends only weakly on when C ^ 4. 
Panel (b) shows how the spectrum depends on the power-law 
index s. Indeed, from equations (1211 1 and (l23T l. it can be seen 
immediately that the nonthermal tail has a broken power-law 
spectrum with indexes (s- l)/2 and .s/2. Panel (c) shows the 
effect of changing 7b. 

The second row shows the dependence of the spectrum on 
Ue, the temperature T (note that we assume that the electron 
temperature is equal to the ion temperature), and the inclina- 
tion angle of the disk /. As the density is increased, the disk 
becomes optically thin at higher frequencies [see panel (d)]. 
The optically-thick portion of the spectrum is blackbody-like 
and is proportional to the projected source size and the tem- 
perature. Increasing the temperature will increase the black- 
body flux linearly. The emission efficiency in the optically 
thin portion is proportional to the square of the tempera- 
ture, which explains the change in shape of the spectrum at 
T — 3 X IO^'K [see panel (e)]. Panel (f) shows the depen- 
dence on the inclination angle / of the disk. On one hand, 
a larger value of / results in a smaller projected source size 
and a small angle between the magnetic field and the line of 
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Fig. 3. — Model fit to the observed broadband spectrum of Sgr A*. The 
horizontal bars show only the extrema obtained from different observations, 
while the arrows at NIR and 7-ray wavelengths give upper bounds on the flux 
densities in the quiescent state. We use T = lO" K and iig = 5x 10* cm^ at the 
outer boundary of the disk. The soHd line is the spectrum computed from our 
emission model with C = 1> * = 3.5, 7b = "o, and cos; = 1/2 from a snapshot 
at the ? = 13 hr of the quiescent simulation. The dashed line coiresponds 
to the perturbed simulation. To demonstrate that this is not a unique fit, we 
compute the quiescent spectrum with a temperature that is 10 times higher 
throughout the disk and f = 10, .s = 2.5 and 71, = 10^ (the other parameters 
are the same as the solid line); the result is shown by the dotted line. 

sight at the blue-shifted side of the disk, implying less efficient 
synchrotron radiation. On the other hand, a larger value of / 
enhances the effect of Doppler boosting. We found that the 
highest peak frequency is reached when cos / =1/2, which, 
obviously, depends on the azimuthal velocity we adopt. 

3. OBSERVATIONAL CONSTRAINTS 

Figure [3] sh ows the observed flux of Sgr A* at different fre- 
quencies (see iLiu et al. 2004, for the instruments and refer- 
ences). The horizontal bars show the maximum and minimum 
fluxes observed so far at the corresponding frequencies. The 
arrows at NIR and 7-ray wavelengths give the upper bounds 
for the quiescent state emission. Fitting these data, one ob- 
tains a (non-unique) set of parameters for the emission model 
and physical conditions of the accretion disk. 

The mm/submm spectrum of Sgr A* has been modeled 
in different ways by various groups. iNarayan et al.l (Il995h 
used thermal emission from an ad vection-dom inated accre- 
tion flow. Mahadevan (1998) and lOzel et"an dlOOO) incor- 
porated the effects of a hybrid electron distribution in the 
same model. .Melia et al. (2001b) considered thermal syn- 
chrotron emission fr om a compact accretion disk, whereas 
iMarkoffetaf] (1200 ll) used the emission of a thermal jet. For a 
magnetic field energy density below equipartition with the en- 
ergy de nsity of the plasma as suggested by MHD simulations 
dHawlev 2000). these studies generally point to an emission 
region of a few Schwarzschild radii with a plasma tempera- 
ture of ^ lO'^K and a density of ~ lO^cm"^. 

Given the low lumi nosity of Sg r A*, the accretion flow is 
radiatively inefficient ( Melia 1992, 1994; iNarayan et al." 19951: 
lYuan et al]r2003, 2004) . The temperature of the plasma there- 
fore should be close to its virial value. As mentioned above. 



we assume for simplicity that the electron and proton temper- 
atures are equal. A realistic treatment of the thermal coupling 
between the protons and the electrons is not feasible in our 
time-dependent simulations. For a given outer boundary ra- 
dius of the accretion disk, the scale height of the disk can be 
estimated from the virial temperature, and we have only the 
gas density left as a free parameter 

We take the outer boundary of the accretion disk to be 
21.5rs in our simulations. We can therefore use the broad- 
band spectrum of Sgr A* to constrain the density at the outer 
boundary, which determines the density normalization of the 
whole disk, the inclination angle of the disk, and other param- 
eters describing the nonthermal population of the electrons 
(specifically, (, s, and 7b). 

The solid line in Figure[3]is an emission spectrum computed 
from a late time snapshot (at f = 13 hr) of our quiescent- 
state simulation (see § |5] for details) with T = lO'^K and 
lie = lO^cm"^ at the outer boundary. The disk inclination an- 
gle is cos; = 1/2, while the emission model parameters are 
( — I, s — 3.5, and 7b — °°. We note that the spectrum peaks 
at 1.2 X lO'^ Hz with a peak flux density of ^ 3 Jy. Ob- 
servations suggest that Sgr A*'s spectrum peaks at 220 GHz, 
slightly above the model predicted value, with a flux density 
of ~ 3 Jy (Marrone et al. 2006). In order to better under- 
stand this difference, we sampled the parameter space to im- 
prove the quality of the fitting. For the chosen snapshot of 
our quiescent-state simulation, we either over-predict the peak 
flux density or have a lower peak frequency as shown in the 
figure. These deviations suggest that the real electron tem- 
perature might be lower than the plasma temperature, as also 
expected on theoretical grounds, while the electron density 
might be higher than in our calculations; this would shift the 
simulated emission peak frequency and flux density closer to 
the observed values. However, quantitative fitting to the spec- 
trum is beyond the scope of this paper, which focuses on the 
variability of the NIR and X-ray emission. We will adopt the 
above model parameters in the following sections. 

Note that, with 7b — °°, we only need two parameters to 
describe the spectrum of nonthermal electrons. Here we ne- 
glect other radiation processes that may dominate in the X- 
rays (i.e., bremsstrahlung and comptonization). Moreover, at 
least part of the quiescent-state X-ray emission from Sgr A* 
is produced near t he capture radius of the acc retion flow by 
a thermal plasma (iBaganoff et al.ll200lL l2003h . As a result, 
the observed X-ray flux needs to be treated as an upper limit 
for the X-ray emission from our model of the small accre- 
tion disk. The above model fitting therefore suggests that for 
C — I, the power-law index can not be less than 3.5. For large 
values of one may adopt a small s without violating the 
X-ray upper limit. These models, however, will predict very 
weak NIR emission. The parameter 7b needs to be included to 
bring the NIR flux close to the observed values during flares. 

The dashed line in Figure |3] corresponds to the perturbed 
simulation discussed in ^ which predicts an X-ray flare with 
a soft spectrum and a flux density 8 times higher than the qui- 
escent level. Further analysis of small flares observed by the 
Chandra and XMM-Newton telescopes are needed to test this 
prediction. As we pointed out above, the values of model pa- 
rameters presented earlier do not constitute a unique fit to the 
spectrum. The dotted line in the figure is another possible fit 
to the spectrum^. The corresponding model parameters are 

* The spectrum is computed by increasing the temperature by a factor of 
ten in our simulation. If we had used Newtonian gravity in our calculation, we 
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Fig. 4. — Gray-scale plots of different quantities in the unperturbed simulation at late time (f = 13 hr). From left to right the panels show the density p, kinetic 
energy density pr jl, magnetic energy density B-/Sit, and the emissivity in the infrared (uq = 1.8 x lO'^^Hz) band with relativistic corrections (cos! = 1/2 in 
this case, see text for details). The polar plots in the top row are vertically-averaged results, while the rectangular plots along the bottom show the variation with 
height in the disk. The asymmetry in the effective emission coefficient comes from the Doppler effect. 



( = 10, s = 2.5, 7b = 10^ and T ^ lO^^K. The values of the 
other parameters are the same as those used to produce the 
soHd Hne. 

4. PHYSICAL SETUP 

We take the mas s of Sgr A* to b e M = 3.4 x lO'^M© 
jS chodel et all 120031: iGh ez et al.| [20031) . From the physical 
conditions in the disk estimated in the previous section, the 
density and temperature at the outer boundary (r 21.5 rs) 
are fixed at po « 1.7 x 10"'^g cm"^ and Tq = 10^ 'K. At this 
temperature, the plasma is fully ionized, so we take the mean 
molecular weight to be /i = 0.5. We also assume the MHD 
equations are valid although the electron distribution has a 
nonthermal component. An absorbing layer is placed between 
20rs < r < 21.5rs to attenuate waves reflecting back into the 
domain. 

In order to focus on the effects of magnetohydrodynamic 
turbulence produced by shearing, we assume a slab geometry 
and neglect the vertical component of gravity. Because the 
disk's mass is negligible compared to that of the black hole, 
we also neglect self-gravity. Inflowing (to the black hole) con- 
ditions are imposed at the inner boundary at r = 1 .5 rs, which 
lies well below the ISCO (« 3rs) in a Schwarzschild geome- 
try. 

At such high temperatures, the electrical conductivity 
within the plasma is high enough that we can neglect any re- 
sistive deviations from ideal MHD (as was done, for example, 
in the global simulations of llgumenshchev & NaravanI [20021 
following the suggestion by Shvartsman 1971). We also ne- 
glect molecular viscosity, which is expected to be insignifi- 
cant compared to the turbulent viscosity arising from Maxwell 
(and Reynolds) stresses in the magnetized plasma. We take 
the initial velocity to be Keplerian in a pseudo-Newtonian 
gravitational potential, v^o = ^ GMrl{r-r%). Shearing leads 
to the development of the MRI (iBalbus & Hawlevl[T99ll) . 

We start with a uniform density and temperature and an 
initial magnetic field B — Bqi. with Bq = 0.3G. A random 
perturbation at the 1 % level is added to the temperature, and 
hence pressure, in order to initiate the MRI. The disk is al- 
lowed to evolve towards a quasi-steady state. Note that a lin- 

could simply rescale the unit of velocity to obtain the answer for a different 
temperature. However, pseudo-Newtonian gravity (or full general relativity) 
sets the characteristic velocity to the speed of light. We cannot rescale the 
temperature without rerunning the simulation. The fit shown here, therefore, 
is only a quantitative demonstration of the uncertainty in our emission model. 



ear mode analysis shows that the modes become stable when 
> 3(ri/vA.z)^, so it is important that the computational do- 
main is large enough to enclose unstable modes for which 
A- > B;{r- r^)^Jlrrj3GMp. For our model, the critical value 
is roughly O.Olr^ at the inner boundary, and \rs at the outer 
boundary. Therefore, we choose the vertical domain to be 
-rs < z < '"s to ensure that the MRI can develop over the 
whole disk. We use 257 x 64 x 32 grid points in our simu- 
lations. The vertical resolution can resolve the most unstable 
wavelength outside the ISCO. The radial and azimuthal grid 
points are chosen so that these two directions are resolved as 
accurately as the z direction. 

5. QUIESCENT STATE 

We carry out the calculation described in ^ for 24 simu- 
lated hours, which corresponds to about 80 orbits at the ISCO. 
The material inside the ISCO is accreted very quickly, in un- 
der an hour. Magnetohydrodynamic turbulence then kicks in 
because of the MRI. The inner disk rotates much more rapidly 
so turbulence is fully developed first in the inner region; the 
MRI in the outer region grows more slowly. Turbulence is 
developed through the whole disk by f = 12 hr (note that the 
period at the outer boundary is around 6 hours), and the disk 
reaches a quasi-steady state thereafter. 

5.1. Quasi-Steady State Solution 

To illustrate the turbulent flow in the disk, we show in Fig- 
ure m gray-scale plots of the density p, kinetic energy den- 
sity pv /2, magnetic energy density B^/Svr, and the NIR (i.e., 
vq= 1.8 X lO^'^Hz) emissivity. The plots show the disk pro- 
file at f = 13 hr of the simulation, after steady state has been 
reached. The NIR emissivity is computed with relativistic cor- 
rections, using the equation 

. /, V ■ ,,,, 

7^0 = 1 -; — ^ T J,^ , (24) 

where we have assumed a 60° inclination angle for the disk. 
The polar plots in the top panels are vertical -averaged quanti- 
ties. The rectangular plots along the bottom show the vertical 
structure at the angles = and tt. Note that the z-axis is 
rescaled to show the vertical structure; the disk itself is much 
thinner than it appears in these plots. 

Because of the high temperature, the density (and the ther- 
mal energy, which is not shown in the plots) is rather smooth. 
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Fig. 5. — The root-mean-square of the different components of the mag- 
netic field. Note that, although we start with a seed magnetic field of 0.3 G 
in the ^-direction, the magnetic field saturates at a value ~ 10 times higher in 
the (dominant) component. 

The region within 3rs is almost empty. Ahhough both the 
density and kinetic energy decrease inside the ISCO, there is 
no significant drop in the magnetic field or the temperature 
there. A weak two-armed spiral pattern, or m = 2 mode, ap- 
pears in the density (and thermal energy). This pattern also 
shows up in the kinetic energy density plot, which is corre- 
lated to that of the magnetic field. The turbulent magnetic 
field shows smaller structure close to the inner disk and larger 
eddies in the out er region. This agrees wi th the standard sta- 
bility criterion of lBalbus & HawlevI (1199 Ih . 

The strong asymmetry of the emission coefficient comes 
from the Doppler shift. Because the plasma moves in a 
counter-clockwise direction (in the polar plots), the emission 
from the left side is blue-shifted. The hybrid emission spec- 
trum decreases for frequencies ly > 10"Hz. Hence, blue- 
shifting raises the emission spectrum and makes the left side 
of the disk much brighter. In the rectangular plots, the patchy 
structure of the emission coefficient for r > 4rs is well corre- 
lated with the magnetic field because j^, oc B^. For the region 
inside 4rs, the low density and low thermal energy, together 
with the gravitational redshift, reduce the emission coefficient 
significantly. 

5.2. Saturated Turbulent Flow 

In Figure |5] we compute the root-mean-square of different 
components of the magnetic field 



jBjrd(f>dz 
J r d<j) dz 



1/2 



(25) 



Although our initial magnetic field has only a vertical com- 
ponent with a seed value of 0.3 G, the dominant compo- 
nent in the quasi-steady state is Z?^. This result is in agree- 
ment with earlier global simulations (see Hawlev 2001). The 
magnitude of the magnetic field (which is dominated by the 
0-component) ranges from 2 G near the outer boundary, to 
the peak value of 4 G near 6rs, and decreases to about 
2.5 G within the ISCO. The magnetic field strength is around 
10% of the equipartition value , consistent with the results of 
iHawlev & Kroiig(l200Tl l2002'). 

Defining an average along the 0- and z-directions as 
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Fig. 6. — Average radial velocity (solid) and magnetosonic speed (dashed) 
profiles thi'ough the disk. The intersection of the two lines defines the mag- 
netosonic point r,ns — 2.4rs. The fluid inside this point loses causal contact 
with the outer part of the disk. 

the magnetic flux along the radial direction is then given by 



if), 



J fr dcj) dz 
J rd(j)dz ' 



(26) 



(27) 



This should vanish for all radii because of the divergenceless 
condition of the magnetic field V • B = and our periodic z- 
axis boundary condition. We compute {Br)ij,z from our sim- 
ulation and obtain {B,)^^ < 10"^^ G at all times. Compared 
to the root-mean-square values of the different components 
of the magnetic field, this shows that our algorithm is able to 
preserve V • B = down to machine accuracy. This property 
is a useful test of our simulation. 

The inner edge of the accretion disk should therefore be de- 
fined as the location where the infalling material loses causal 
contact with the rest of the disk. Hence, it may be functionally 
defined to be the magnetosonic point, at which the radial ve- 
locity matches the (fast) magnetosonic speed Cf = •\/cs2 + ca^, 
where Cs is the sound speed and ca~B/ y/Airp is the Alfven 
speed. Below this radius, the material loses causal contact 
with the rest of the disk at larger radii. We plot in Figure |6] 
the root-mean-square of the radial velocity and magnetosonic 
speed as functions of radius. The intersection indicates the lo- 
cation of the magnetosonic point, which is around 2.4rs. The 
Keplerian period at this radius is 10.76 min, in contrast to the 
17.19 min period at the ISCO. 

6. FLARES 

Earlier hydrodynamic and MHD si mulations of the large- 
scale (r > 20rs) accretion flow (e.g. iRuffert & M eHa 1994; 
ICoker & Melia'4997t llgumenshchev & Narayanll2002. ,2003 : 
Cuadra et al. 200lF have indicated that the inflow onto Sgr 
A* is not smooth. The time scales of these simulations range 
from 1 to lO''^ years. Given the low-density conditions in 
the surrounding medium, only parcels of plasma with rel- 
atively sma ll specific angular momentum find their way to 
small radii dCoker &' Melia 1997). The lack of strong accre- 
tion apparently inhibits the formation of a large continuous 
disk extending out to the Bondi-Hoyle capture radius (around 
lO^rs). Instead, the compact disk in Sgr A* appears to be 
accreting clumps of plasma th at "rain" inwards fr om all di- 
rections (Falcke & Melia 1997; Coker. MeUa. & F alcke 199^ 
iTag ger & M elia 2006). 

Each clump circularizes at a radius corresponding to its spe- 
cific angular momentum. Some clumps presumably reach as 
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Fig. 7. — Same as Figure|4]but for the perturbed simulation. 
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far in as the ISCO; most of the others probably merge with 
the disk at larger radii. To model the impact of such a clump 
falling onto our quasi-steady inner disk, we simulate its effect 
on the luminosity of the disk by introducing a density pertur- 
bation in the saturated disk. 

6.1. Density Perturbation 

The perturbation is introduced into the quiescent simulation 
after quasi-steady state is reached, at f = 12 hr. The perturba- 
tion is Gaussian in density 



bp = 2poexp 



r^ + rQ + 2rroCOS( 
la 



(28) 



where po = lO^mn cm"-', ro = 6rs and a = 2r%. We assume 
that the perturbation has a temperature of only lO^^K, much 
cooler than the disk, to avoid a strong increase in pressure. 
This allows the clump to move with the plasma in the disk in- 
stead of propagating out as sound waves. The internal energy, 
therefore, is raised by 



bE ^ Sp- 



3k. 



X 10^°K. 



(29) 

For simplicity, we also assume the clump has zero momen- 
tum. Hence the velocity is slowed down by 

-Sp 



6V: 



(30) 



p+Sp 

The simulation is carried out for 10 hours, until the accre- 
tion rate becomes comparable to that of the quiescent state 
again. In Figure|7] we plot the density, kinetic energy density, 
magnetic energy density, and the effective NIR emission coef- 
ficient at f = 13 hr. The various physical quantities are shown 
at the same time and with the same gray-scale as in Figure |4] 
for comparison. The perturbation is sheared out and forms 
a one-armed spiral pattern in the disk, which also raises the 
kinetic energy. The correlation between the magnetic energy 
density and the gas density is less clear in the plot. The per- 
turbation seems to produce a few strong magnetic spots, but 
it also weakens the m = 2 mode. There is a prominent feature 
in the effective emission coefficient; this "hot spot" arises due 
to both the density and the temperature enhancements. 

6.2. Accretion Rate 

Before looking for quasi-periodic signals in the lightcurve, 
we first study variations in the accretion rate. In Figure |8] we 



plot the accretion rate of our simulation at r = 3rs as a func- 
tion of time. The solid line corresponds to the unperturbed 
simulation, and the dashed line corresponds to the accretion 
rate after we introduce the density perturbation. The strong 
peak in the first hour is due to our initial condition of a uni- 
form disk. The material inside the ISCO falls into the black 
hole immediately after the simulation begins. Although the 
accretion rate never settles to a constant value, it converges 
to a value around 2 x lO'^g s~', consistent with the value in- 
ferred from earlier semi-analytic treatments. The weaker peak 
appearing right after t = 12 hr is caused by the density pertur- 
bation. 

Although we do not show it in the plot, the accretion rate 
M at radii r > 6rs oscillates between positive and negative 
values (it actually ranges from -6 x lO'^ to 1 x lO'^g s~'). 
These oscillations in M have well-defined periods of about 1 
hour, though they are not related to the orbital period. They 
are simply pressure-driven. 

6.3. Disk Images 

Assuming an inclination angle of 60° and using the param- 
eters obtained in ^ we integrate the transfer equation along 
parallel rays and obtain ray-traced images of the accretion 
disk around Sgr A*. Compared to the emissivity plots in Fig- 
ures |4]and|7] these images represent the observation more di- 
rectly because they are numerical solutions of the radiative 
transfer equation. They also take into account the projection 
effects of the disk. 

In Figures |9] and [TO] we present the ray-traced im- 
ages of our unperturbed and perturbed simulations, respec- 
tiv ely. Note, however, that u nlike the images provided 
in iFalcke. JVIelia. & Agoll (|2000|) and iBromlev. Melia. & Liul 
(I2001h . in which the effects of interstellar scattering and the 
diffraction-limited finite resolution of the telescope array were 
taken into account to produce realistic images observed at 
Earth, these are meant only to show the emission character- 
istics at the source. 

The images show the disk at mm (10^' Hz, left), NIR 
(^Qi4.25 middle), and X-ray (lO'** Hz, right) wavelengths, 
from f = 12 hr (top) to 14 hr (bottom). The K-band is cen- 
tered on 2.2 microns and has a bandwidth of 0.48 microns. 
We approximate it by the frequency 1.78 x lO'^Hz, which is 
equal to 10^'^'^^Hz. The color scales are chosen to enhance 
the contrast in the emission and are fixed for each frequency. 
The value of /Vq given in the plot is the flux density observed 
at Earth from the whole image, according to equation (|2]). 
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Fig. 8. — The variation of the accretion rate M at r = 3rs with time. The solid line corresponds to the unperturbed simulation. The peak that appears in the 
first hour is due to our initial condition of a uniform disk. Matter below the ISCO falls into the black hole very quickly after the simulation begins. The dashed 
line after t = 12hr is the result from the simulation with a density perturbation (see The accretion rate rises for about half an hour, then falls back to the 
quasi-steady level. 



From Figure |9] we already see that the accretion disk is 
highly variable. The flux densities of the unperturbed disk at 
12 hr 20 min and 13 hr 40 min are different by a factor of 
two in both the NIR and X-ray bands. They are only differ- 
ent by about 70% in the mm/sub-mm bump. These variations, 
due to the bright structure at lOrs [panel (bO)], are consistent 
with the observed frequency dependence of the flux variabil- 
ity. This suggests that the turbulent structure itself can only 
generate factor-of-two fluctuations. However, although it is 
not likely to happen, because we have only run the simulation 
for 24 hours, we cannot rule out the possibility that the turbu- 
lence can generate flares that increase the disk brightness by 
a larger factor, though much more rarely than these smaller 
variations. 

The first row of Figure [TO] for the perturbed state is iden- 
tical to that of Figure |9] because the perturbation is added at 
= and t — 12 hr The flow shears out the perturbation and 
forms a one-armed spiral pattern as shown in Figure |71 Not 
surprisingly, there is a spiral pattern in the corresponding im- 
ages. At f = 12 hr 20min, we start to see a hot spot appear 
near the inner edge of the disk. This hot spot is much brighter 
than the bright structure seen in the unperturbed simulation. 
The hot spot keeps moving outward due to the leading spiral 
pattern of the perturbation. One hour after the perturbation 
is introduced, the disk reaches its brightest moment when its 
flux is 20% higher than its quiescent value in the radio and 
500% brighter in the NIR and X-ray bands. 

In addition to the bright one-armed spiral, there are other 
bright features in the inner part of the perturbed disk. For ex- 
ample, there are two hot spots in the f = 13 hr image. The 
inner hot spot moves much faster than the spiral pattern, and 
appears in the image about every half hour after the pertur- 
bation is introduced. These corotating features suggest that 
QPOs may be found in these simulations. 

6.4. Light Curves and QPOs 

In Figure [TTl we show the NIR (upper) and X-ray (lower) 
lightcurves obtained by integrating the flux densities over the 
disk images. The solid and dashed lines correspond to the 
quiescent and flare state simulations, respectively. 

There are several important differences between the 
lightcurves shown in Figure [TT] and the accretion rate plotted 
in Figure[8] First, note that the initial peak in the accretion rate 
in Figure 10 is a transient phenomenon that occurs because the 
density relaxes from an initial constant value to a steady-state 
solution. At the same time, the MRJ has not yet amplified 
the magnetic field to produce significant synchrotron radia- 
tion. Second, the accretion rate is noisier than the lightcurves 



because it is a surface integral at some specific radius (at 3rs 
in our case), while the flux density is a volume integral over 
the whole computational domain. Features with a width A in 
the accretion disk pass through a specific radius over a time 
scale A /v/, however, they can affect the lightcurve for a much 
longer time (rmax - '•min) Ar- 

The flare near f = 12 hr in the perturbed lightcurves lasts 
much longer than the corresponding peak in the accretion rate. 
The NIR (and X-ray) flux density reaches its peak att = 13 hr, 
when the accretion rate at 3.0rs has dropped back to the qui- 
escent level. These results show that the flux densities are 
not necessarily correlated with the accretion rate. Indeed, our 
simulations indicate that the perturbation not only raises the 
density trivially, but is also able to change the structure of the 
turbulence. This kind of distortion takes about an hour to de- 
velop, and merges with the original turbulent disk within a 
second hour 

Comparing the lightcurves to the disk images, we can iden- 
tify the second peak of the flare with the hot spots shown 
in panels (b3) and (c3) in Figure [10] It is directly related 
to the spiral structure in the disk (see Figure [T])- However, 
the oscillations during the flare, which persist all the way to 
about f = 20 hr, come from other features developed by the 
perturbation — for example, the smaller hot spot located at the 
inner boundary in panel (b3) of Figure [TOl 

The evolution of these small fluctuations suggests that 
there should be a periodic (or quasiperiodic) signal in the 
lightcurves. In Figure [121 we over-plot the power spectral 
densities (PSD; multiplied by the frequency) in the NIR band 
in both the quiescent and flare state simulations. The curves 
are normalized at higher frequencies for easy comparison of 
the two. The PSDs of the X-ray lightcurves are almost iden- 
tical to (and hence overlapping with) the NIR PSDs because 
both of them are produced by the nonthermal electrons. The 
two arrows indicate the frequencies associated with the orbital 
period at the ISCO and the magnetosonic point (at 2.4rs). 

Although there is no significant period in the quiescent 
state, a QPO-like feature is found in the flare simulation. The 
quasiperiod is around 1 1 min, essentially the Keplerian pe- 
riod at the magnetosonic point. This is, to our knowledge, the 
first QPO seen in magnetohydrodynamic turbulent disks. If 
such a density perturbation (due to "raining" plasma) is in- 
deed the cause of such flares in nature, the QPOs observed 
from Sgr A* may be used to set a lower limit on the orbital 
period at ISCO, and hence provide an estimate of the spin of 
the Galactic center black hole. This idea can be applied to 
other low-luminosity AGN as well. 

7. LIMITATIONS 
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Fig. 9. — Ray-ti'aced images of the unpeiturbed simulation of tlie accretion dislc in Sgr A*. The plots show the images at the thermal peak (left), infrared 
(middle), and X-ray (right) bands from I = 12 hr (top row) to 14 hr (bottom). The color scale for each frequency is fixed so that the reader can see the variation 
of the image with time and compare with Figure flOl for the perturbed simulation. Here, is the flux density that would be observed at Earth (from the whole 
image). Note, however, that these images do not include the distortions produced by interstellar scattering and the finite telescope resolution due to diffraction, 
and are meant only to convey a sense of the emission geometry at the source. 
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Fig. 10. — Same as Figure|9]but for the perturbed simulation. 
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Fig. 11. — NIR and X-ray lightcurves. The solid lines are computed from the quiescent state, and the dashed lines are from the flare simulation. Note 
that the peaks appearing at the 13th hour ai'e much broader than the peak seen in the accretion rate. This indicates that the perturbation actually disturbs the 
magnetohydrodynamic disk, instead of trivially raising the density. There are some oscillations on top of the peak, which persist all the way to the 20th hour 
These oscillations originating from hot spots in the disk produce a QPO in the lightcurve. 



Although for the first time a QPO is seen in magneto- 
hydrodynamic simulations of accretion disks without invok- 
ing special conditions such as large-scale magnetic fields, 
the QPO period is shorter than that observed from Sgr A* , 
which is usual l y long er than 17 minutes dGenzel et alj|2003l : 
iBelangeretaf] 120061) . This may be partially due to the 
pseudo-Newtonian potential we have adopted here. In a 
Schwarzschild geometry, the orbital period at the ISCO for 
a 3.4 X lO^Mp black hole is 25.7 minutes. If the magne- 
tosonic point is still at 2Ars in such a potential, the cor- 
responding period is 18.4 minutes, which is slightly longer 
than the period observed in the NIR band but compatible 
with the X-ray observations ("Belanger et al.l 12006 ). A more 
complete sirnulation incorporating general relativistic effects 
(e.g., iHirose et al.il2004j) and relati vistic ray tracing calcula- 
tions ( iBromlev. Melia. & Liul20()ll) is needed to support more 
quantitative comparisons with observations. 

We also assume a periodic boundary condition along the 
vertical direction and ignore the vertical component of grav- 
ity. Therefore our simulations are only 2. 5 -dimensional, 
in the sense that they are able to correctly simulate three- 
dimensional turbulence but cannot capture many other three- 
dimensional features. For example, the geometry does not 
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Fig. 12. — The PSD (multiplied by the frequency) of the lightcurves. The 
curves are normalized so that the flat noise spectra at high temporal frequen- 
cies match each other There is no significant period during the quiescent 
state. In the perturbed simulation, a QPO appears at 11 min, which is the 
orbital peiiod at 2.4rs, the magnetosonic point. 



allow an outflow along the vertical direction. Hence, it cannot 
be used to determine whether there are jets associated with 
the hot accretion disk, which may contribute to nonthermal 
radio emission from Sgr A* (Liu & Melia 2001) and explain 
the ob served cor relation between X-ray flares and radio out- 
bursts dZhao et al.. 2004) . We also cannot see magnetic-field- 
dominat ed funnels as revea led in global relativistic MHD sim- 
ulations (^Hirose et al.l l2004ft . These funnels may be responsi- 
ble for the acceleration of energetic protons near the black 
hole ( Liu et alj |2006b) and can play an important role in driv- 
ing Poynting-flux-dominated outflows, which are proposed to 
power large-scale jets in black hole accretion systems. 

The saturation level of our ma gnetic field agre es with 
the glob al simulation of |Hawley & Krolik ( 1200 IL [2002). 
iHawlev et al.l (119951) pointed out that the saturation level of 
the magnetic field in shearing box simulations depends on the 
box size. The verti cal size of our sim ulation domain therefore 
is a free parameter ' Stone et al] (1199 6) also showed that there 
are more complicated structures in the vertical direction when 
the effects of gravity are included. A fully three-dimensional 
simulation is needed to resolve these issues. 

To compare observations with MHD simulations, one also 
has to convert the simulated disk characteristics into a radia- 
tion flux spectrum. Observations have already provided good 
constraints on the emission processes. However, it is not clear 
how electrons are energized in the turbulent plasma, which is 
related to the fundamental problem of non-linear turbulence 
dissipation and electron-ion coupling. In the case of rela- 
tivistic collisionless plasmas considered here, we lack even 
the appropriate tools to handle this subject self-consistently. 
There have been only some phenomenologic al models devel- 
oped to address this p roblem quantitatively dLiu et al.ll2004 
l2006al:lBittner et al.l2006.) . In this paper, we have adopted one 
of the simplest emission models appropriate for the accretion 
flow in Sgr A* and have used observations to constrain the 
model parameters. Although the QPO signal is not very sen- 
sitive to the model parameters as long as there is a higher en- 
ergy electron population above the thermal background, one 
should exercise caution with more quantitative comparisons. 

Besides the points mentioned above, there are additional 
physical processes we can include to improve our sim- 
ulations and quantitative results. For example, thermal 
bremsstrahlung, and inverse Compton scattering of the syn- 
chrotron photons by the hot electrons may provide a signifi- 
cant contribution to the X-ray emission (lNarayanetalJll995h 
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especially during big flares dLiu et al.ll2004l) . 

8. CONCLUSIONS 

In fliis paper we have simulated a nonradiative MHD accre- 
tion disk around Sgr A* using a pseudo-spectral algorithm. 
The results are consistent with earlier th ree-dimensional 
MHD simulations of the MRl ( iHawlev & Bal bus 2002), even 
though we have assumed a different geometry. Our simula- 
tions reached a quasi-steady state after the MRl-driven turbu- 
lence has fully developed. 

We computed the long-wavelength spectrum from our sim- 
ulation by using a hybrid emission scheme with a thermal 
background of electrons plus a high-energy nonthermal bro- 
ken power-law tail. We then used the broadband quiescent- 
state spectrum to calibrate the model parameters. In the quasi- 
steady state we found that the source is variable with an am- 
plitude in flux density that increases with photon frequency, in 
agreement with observations. However, the flux density var- 
ied by only a factor two, failing to account for the observed 
NIR and X-ray flares from Sgr A*. 

Motivated by earlier st udies of the large-scale accretion pro- 
cesses in Sgr A* (e.g., ICoker & Melial 119971; ICuadra et all 
12001 . we introduced a density perturbation to the quasi- 
steady state disk to see whether they could produce even 



bigger flux density variations. We found that bigger flares 
were produced in both the NIR and X-ray bands, as well 
as a QPO in the lightcurves in association with the magne- 
tosonic point below the ISCO. The fact that the simulated 
QPO is associated with the magnetosonic point below the 
ISCO also requires a reexamination of previous estimates of 
the black hole spin based on fl are observations (iGenzel et alj 
l200llAschenbach et al.ll2004 . 
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